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Spin-dependent electronic transport properties of Fe-terminated zig-zag graphene nanoribbons (zGNR) have 
been studied using first-principles transport simulations. The spin configuration of proposed zGNR junction 
can be controlled with external magnetic field, and the tunneling junction show MR>1000 at small bias and 
is a perfect spin- filter by applying uniform external magnetic filed at small bias. 

PACS numbers: 72.80.Vp,73.63.Rt,73.22.Pr 
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Graphene, a single layer of graphite honeycomb lat- 
tice, has attracted great attention since its discovery 
since 2004^. The material exhibits peculiar electronic 
structure and large mean free path, and has been there- 
fore considered a potential replacement for silicon in the 
future electronics^. Son et al. studied the electronic 
structure of graphene nanoribbons (GNR), and proposed 
that the zig-zag GNR (zGNR) serve as electrically con- 
trollable spin-valve^. This result quickly draws atten- 
tion and the study of GNRs becomes the focus of the 
graphene research. Various methods to prepare GNRs 
or to enhance their qualities have been proposed^ 
the edge states and stabilities have been explored^^-; 
the doping and edge-termination properties have been 
studied^' ^-iii^ii2r,^^; and possible applications in electronic 
devices have been discusse d-'^^i^^" — . In this letter, we 
show a possible implementation of GNR based spintron- 
ics device, which shows large magnetoresistance (MR) as 
well as perfect spin- filtering effect. 

The proposed tunneling junction is modelled with two 
semi-infinite Fe-terminated zGNRs with hydrogen ter- 
mination at the ends (Fig. [D^a)). Previous studies 
have shown that the transition metal terminated zGNRs 
(TM-zGNR) exhibits both antiferromagnetic (AFM) and 
ferromagnetic (FM) long-range ordering like the pres- 
tine zGNRs and the hydrogen-terminated zGNRs (H- 
zGNR)^^. Therefore, at least three distinct spin states 
can be realized with the proposed junction: 1. edges 
within each lead are ferromagnetically coupled, while 
edges across the leads are antiferromagnetically coupled 
(Fig. mb)); 2. edges within each lead are antiferromag- 
netically coupled, while edges across the leads are ferro- 
magnetically coupled (Fig. [I{c)); 3. edges within each 
leads as well as across the leads are ferromagnetically 
coupled (Fig. md)). In the following text, these states 
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FIG. 1. (color online) Geometry (a) and spin configurations 
(b-d) of the model system. The dashed line enclosed region is 
the scattering region in simulations. The arrows on the atoms 
indicates the direction of local magnetic moment. 



will be referred to as the AP state, the AF state, and the 
PR state, respectively. Experimentally, the AF state is 
the ground state without any external B-field, while an 
uniform external B-field will turn the system into the PR 
state. The AP state could be induced if external B-fields 
in opposite direction are applied to the two leads. 

In order to model the transport properties of the pro- 
posed tunneling junction at different states, we employed 
non-equilibrium Green's function (NEGF) method com- 
bined with density functional theory (DFT). The struc- 
ture of the modelled system was fully relaxed using the 
plane- wave basis DFT code PWSCFi^a, with the ultrasoft 
pseudopotential method. The energy cutoffs for plane- 
wave and density were chosen to be 40 Ry and 320 Ry, 
respectively. The transport properties was then obtained 
with the optimized geometry using the SIESTA code^. 
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FIG. 2. (color online) Transmission coefficient G(E) of the 
tunneling junction at AP (upper panel), AF (middle panel), 
and PR (lower panel) states at zero-bias. Please be noted that 
the conductance scale of the upper panel is two magnitudes 
larger than the other two. The Fermi level is shifted to (the 
green dashed line). 

with the norm-conserving pseudopotential method and 
atomic orbitals. The real-space mesh grid size was chosen 
to represent 200 Ry equivalent plane-wave energy cutoff, 
and single- C basis with optimized parameters were em- 
ployed to avoid the singularity problems in calculating 
the lead Green's function. In all calculations, Perdew, 
Berke, and Ernzerhoff parameterization (PBE) to general 
gradient approximation (GGA)^^ was used to calculate 
the electron exchange-correlation. 

We first analyze the transport properties of this GNR 
junction. We show the transmission coefficient G{E) of 
the GNR junction at different states in Fig. [2l Due to 
the symmetry of the junction and the spacial spin con- 
figuration, the a-spin and /3-spin channels are degenerate 
for both AP and AF states, as shown in the figure. For 
the tunneling junction at the PR or AF states, since both 
leads are exactly the same, the electron band energies of 
both leads in either spin channel perfectly match each 
other, and therefore the transmission coefficient G{E) of 
these states are two magnitudes larger than G{E) of the 
AP state. As the state of the tunneling junction could 
be controlled with external magnetic field, it is therefore 
expected that the tunneling junction will show significant 
magnetoresistivity behavior. For the tunneling junction 
at the PR state, it can further be noticed that the degen- 
eracy over spin channels is lifted, and a large conductance 
peak shows up at the Fermi level Ep for the /3-spin chan- 
nel, while the a-spin channel does not open until around 
0.4 eV. Thus, it is expected that the junction could serve 
as a perfect spin-filter at PR state. 

To further illustrate the possible MR and spin-filtering 
behavior, we employed the following equation 

1 = j ^-h^^E) {fUE)-fn{E))dE 




0.5 1 1.5 

Bias (V) 

(a) 




I ^ ^ ^ 1 

0.5 1 1.5 

Bias (V) 

(b) 

FIG. 3. a) Magnetoresistance of PR/AP and AF/AP states; 
and b) current-polarization at PR state. The y-d^xis in panel 
a) is in logarithm scale. 

to calculate the current at small bias. The obtained cur- 
rent is then used to calculate the MR between PR/AP 
and AF/AP states (FIG. [3fa)). As expected, the junc- 
tion shows MR order of 10^ between either PR/AP states 
or AF/AP states for a wide range of bias voltage. In the 
AF/AP case, the small conductance hump around Ef 
at the AP state causes the MR initially smaller than 1. 
However, as the conductance peak of AF state is close to 
Ef and its scale is two magnitudes larger, the current at 
AP state quickly become dominant and the MR reaches 
over 1000 from 0.3 to 1 V. For the junction at the PR 
state, the obtained spin-current is also used to calculate 
the current polarization ^ defined with 

^ ^ 1(3 - I a 

where a and f3 denotes majority and minority spin, re- 
spectively. As shown in FIG.[3]^b), the spin polarization 
of the current is almost 100% for small bias (F < 0.5 V). 

In order to understand the physics behind the GMR 
and spin- filtering behavior of the Fe-terminated GNR, we 
have performed density of states (DOS) calculations. In 
either AP or AF states, both spin channels of the scat- 
tering region are degenerate as shown in FIG. lU^a) and 
mb) due to the spacial symmetry, in consistency with the 
transmission coefficient data. For the tunneling junc- 
tion at AP state, although the scattering region DOS 
show large peaks around Ef for both spin channels, the 
left/right lead DOS are almost negligible for a-//3-spin 
respectively (the figure shows only the left lead DOS, the 
right lead DOS is the same except that the spin channels 
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FIG. 4. Density of states of left lead and scattering region at a) AP, b) AF and c) PR states. The upper panels are for a-spin; 
while the lower panels are for /3-spin. 



needs to be swapped). Furthermore, the spin transfer for 
this state is also forbidden due to the mismatch between 
orbital symmetries at both sides, similar to the case in the 
pure zGNR device^^. Therefore, neither spin can easily 
tunnel through the junction, and thus the transmission 
coefficient at AP state is orders of magnitudes lower than 
other two. The scattering region DOS at AF state show 
large peaks around 0.1 eV, where small humps are also 
present in the leads DOS. This feature implies resonant 
tunneling through the junction, which leads to the large 
transmission coefficient peak at AF state around 0.1 eV 
(middle panel of FIG. [2j). Finally, for the tunneling junc- 
tion at the PR state, both leads and the scattering region 
have the same orbital symmetry as well as aligned energy 
levels^^, thus a perfect transmission is expected. How- 
ever, as both leads and scattering region show negligi- 
ble DOS for the a-spin channelU^c), therefore the a-spin 
conductance is greatly suppressed. In contrast, for the 
/3-spin channel, the scattering region DOS show a large 
peak at E^^ where a small DOS hump is also present for 
the leads DOS. Thus, resonance tunneling takes place 
at small bias for the /3-spin. Hence, the junction show 
significant spin- filtering effect at small bias. 

In conclusion, we have performed ffist-principles trans- 
port simulations on Fe-terminated GNR junctions, whose 
spin configurations may be controlled by applying exter- 
nal magnetic field. The junction show MR over 1000 at 
small bias for a large bias range, and is a perfect spin- 
filter at PR state at small bias. The behavior of the junc- 
tion is controlled by the DOS of both leads and scattering 
region as well as the associated orbital/band symmetries. 
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